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Abstract 

We use our recently developed functional bosonization approach to bosonize 
interacting fermions in arbitrary dimension d beyond the Gaussian approxi- 
mation. Even in d = 1 the finite curvature of the energy dispersion at the 
Fermi surface gives rise to interactions between the bosons. In higher dimen- 
sions scattering processes describing momentum transfer between different 
patches on the Fermi surface (around-the-corner processes) are an additional 
source for corrections to the Gaussian approximation. We derive an explicit 
expression for the leading correction to the bosonized Hamiltonian and the 
irreducible self-energy of the bosonic propagator that takes the finite curva- 
ture as well as around-the-corner processes into account. In the special case 
that around-the-corner scattering is negligible, we show that the self-energy 
correction to the Gaussian propagator is negligible if the dimensionless quan- 
tities (^) d i ? o[l + ^ ? o]~ 1 ^|^"| are small compared with unity for all patches 
a. Here q c is the cutoff of the interaction in wave-vector space, kp is the 
Fermi wave- vector, fi is the chemical potential, Fq is the usual dimensionless 
Landau interaction-parameter, and v a is the local density of states associated 
with patch a. We also show that the well known cancellation between vertex- 
and self-energy corrections in one-dimensional systems, which is responsible 
for the fact that the random-phase approximation for the density-density cor- 
relation function is exact in d = 1, exists also in d > 1, provided (1) the 
interaction cutoff q c is small compared with Jzf, and (2) the energy disper- 
sion is locally linearized at the Fermi surface. Finally, we suggest a new 
systematic method to calculate corrections to the RPA, which is based on the 
perturbative calculation of the irreducible bosonic self-energy arising from the 
non-Gaussian terms of the bosonized Hamiltonian. 

PACS numbers: 05.30Fk, 05.30.Jp, ll.10.Ef, 71.27.+a 
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I. INTRODUCTION 



One of the most powerful methods to analyze systems of interacting fermions in d = 1 
dimension is the bosonization approach jp. Over the past 30 years numerous interesting 
results have been obtained with this essentially non-perturbative method, and the study 
of one-dimensional fermions via bosonization techniques continues to be a field of active 
research. Because interacting fermions in d = 1 are not Fermi-liquids ||, conventional many- 
body perturbation theory is not applicable to these models. The success of bosonization in 
d — 1 hinges on the linearization of the energy dispersion at the two Fermi points. As 
soon as the curvature of the energy dispersion is taken into account, bosonization maps 
an interacting Fermi system onto an effective interacting Bose system, which in general 
cannot be solved exactly. Of course, one can circumvent this problem by chosing the energy 
dispersion of the Fermi system to be linear by definition || (Tomonaga-Luttinger model), 
but then the physical relevance of this model needs further justification, because in realistic 
materials the energy dispersion is never exactly linear. 

In higher dimensions conventional many-body perturbation theory is at least consistent, 
provided the interaction is not too singular Because within the framework of perturba- 
tion theory there is no sign for the breakdown of Fermi-liquid theory, there was no need to 
develop non-perturbative methods for interacting fermions in dimensions higher than one. 
Luther's pioneering attempt to describe higher dimensional Fermi-surface degrees of free- 
dom by bosonization was not further pursued |§. The discovery of the high-temperature 
superconductors and Anderson's suggestion P,|lT)|| that the normal state properties of these 
materials are a manifestation of non-Fermi-liquid behavior in d = 2 has revived the interest 
to develop non-perturbative methods for analyzing interacting Fermi systems in d > 1. Be- 
cause of the success of bosonization in d — 1 it was natural to further develop this approach 
in higher dimensions. Recent progress in this field is based on ideas due to Haldane flIT] , |T2"f , 
who realized that Luther's construction || must be generalized such that the direction of the 
net momentum of particle-hole pairs is not restricted. (In Luther's approach particle-hole 
pairs carry a net momentum that is strictly normal to the local Fermi vector). 



For later reference it is useful to summarize at this point the basic features of Haldane's 
"tomographic" construction. Consider the many-body Hamiltonian 

H = H + H int , (1-1) 

^o = E e k4a k , (i-2) 

k 

Him = JttY, /£ k '4+q4'-qCk'C k , (1.3) 
ZV qkk' 

where Ck annihilates an electron with wave-vector k, V is the volume of the system, and 
/q k are generalized Landau-parameters that depend not only on the momentum transfer, 
but also on the momenta of the incoming particles. For simplicity we have suppressed the 
spin label. The spin is easily taken into account by defining k and k' to be collective labels 
for wave-vector and spin. We are implicitly assuming that only a single band with energy 
dispersion of the underlying lattice model lies in the vicinity of the Fermi energy, and that 
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the effect of the other bands is either completely negligible or can be taken into account via 
the definition of effective parameters e k and /q k . The Fermi surface is the d — 1-dimensional 
hypersurface in k-space which satisfies 

& = e k -p = , (1.4) 

where \i is the chemical potential, measured with respect to the bottom of the conduction 
band. In d = 1 the Fermi surface consists of two distinct points ±kp, where kp is the 
Fermi wave- vector, while in d > 1 it is a d — 1-dimensional manifold, the topology of which 
depends on the form of £ k . The crucial step in Haldane's construction is the subdivision of 
the Fermi surface into disjoint patches with volume A d_1 . The precise shape of the patches 
and the size of A should be chosen such that within a given patch the curvature of the Fermi 
surface can be locally neglected, so that the variations of the local normal vector can be 
ignored within a patch. Note that in this way sufficiently flat Fermi surfaces can be covered 
with only a small number of patches. The extreme case is the Fermi surface of an array of 
chains in d = 3 without inter-chain hopping, which consists of two parallel planes that can 
be identified with the patches. Obviously, in this case two patches are sufficient to cover the 
entire Fermi surface flT3| , |H|| . 

We introduce a label a to enumerate the patches in some convenient ordering and denote 
the patch labelled a by K%. Note that by definition K% is a subset of the Fermi surface, i.e. 
it is also ad — 1-dimensional hypersurface. To describe all degrees of freedom, we extend 
each patch K% into a d- dimensional "squat box" K^ x of radial height A, such that the 
union \J a K% A agrees with the wave- vector space necessary to label all degrees of freedom in 
the system, see Fig.|l]. By introducing a radial cutoff A, we are implicitly assuming that all 
states with wave-vectors outside a shell with thickness A around the Fermi surface have been 
integrated out, and that the associated finite renormalizations have been taken into account 
via the proper definition of the bare parameters e k and / kk in Eqs. |1.2| and |1.3| . Of course, 
this construction is only useful if the physical quantities of interest do not depend on the 
choice of the cutoffs A and A. For sufficiently smooth Fermi surfaces this is indeed the case 
as long as we are interested in correlation functions at separations |r — r'| ^> max{A^ 1 , A -1 }, 
and the interaction is dominated by momentum transfers |q| <C min{A, A} ||16|| . However, 
the slightly vague picture of integrating out the high-energy degrees of freedom without 
really doing this calculation is not very useful if we are interested in the precise numerical 
value of some physical quantity. To describe all degrees of freedom, we should therefore 
remove the radial cutoff A and work with boxes K% that are constructed such that \J a 
covers not only a thin shell around the Fermi surface, but also the high-energy degrees of 
freedom. Of course, if we would like to linearize the energy dispersion around the Fermi 
surface, we should keep the radial cutoff finite, because the linearization is only good close to 
the Fermi surface. However, our bosonization approach does not rely on the linearization, so 
that in some cases we can remove the radial cutoff and can dispose of the picture of having 
integrated out the high-energy degrees of freedom. For example, to discuss the homogeneous 
electron gas in d = 2 we take e k = J^, / kk ' = an d choose the boxes K% shown in Fig.||. 
We shall come back to the cutoff problem in the radial direction in Sec. [lV C) . 

Because the union of all K a agrees by construction with the total relevant k-space, we 
have 
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E© a (k) = l , (1-5) 

a 

where 



For simplicity we have omitted the cutoff label of K a . In the relevant Hilbert space Eq. |1.5| is 
correct for both type of boxes K% and K% A discussed above. Let us denote by k a a vector on 
the Fermi surface (£ k <* = 0) that points to the (suitably defined) center of patch a. The set 
of vectors {k a } define the origins of local coordinate systems on the Fermi surface. Such a 
collection of coordinate systems is also called an atlas. Let us now measure the wave- vector 
k locally with respect to the origin of the appropriate coordinate system. Using Eq.[L5], we 
may write 

£ k = £e«(k)£ k = £e«(k)^_ kQ , (i.7) 

a a 

where the quantity £° is defined by 

The crucial step in the conventional bosonization approach is the local linearization of the 
energy dispersion within a given patch, 

q = v a -q + 0(q 2 ) , (1.9) 

where 

v a = V k e k | k=kQ (1.10) 

is the local Fermi velocity at patch a, and we measure frequencies in units of energy, which 
amounts to formally setting % = 1. Note that such a linearization would not be possible 
in a rigid coordinate system. It is important to stress that this linearization should be 
made around the true chemical potential of the interacting many-body system [|TjJ, and 
that v a should be identified with the renormalized local Fermi velocity, taking effective 



mass renormalizations into account ||16|| . 



The operator-approach to bosonization follows then closely the usual construction for 
one-dimensional systems M. This approach has been further developed by Houghton, Kwon, 



Marston and Shankar []17|- [pOjl , and independently by Castro Neto and Fradkin [glj]- |23 
The fundamental objects are the local density operators associated with the patches, 



P. 



« = ]Te a (k)5 k £ k+q . (1.11) 



In terms of these operators the interaction-part of H can be written as 

^=^EE/r' : «' : > ( Li2 ) 



1 aa' 
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where : . . . : denotes normal ordering, and it is assumed that the variations of /^ k ' are 



q 



negligible if k and k' are restricted to given patches, so that it is allowed to introduce the 
coarse-grained interaction function 

w _ E kk ^(k)9«'(kQ/f ' 

q Ekk'0 a (k)e«'(k') ■ 1 ' 1 

By considering the commutation relations of the operators with each other and with H 
in the restricted Hilbert space of states with wave-vectors close to the Fermi surface, it is 
possible to show that the approximately obey bosonic commutation relations, and that H 



is approximately quadratic in the local density operators, so that it is exactly solvable |12 
0- li- 



lt should be mentioned that an alternative non-perturbative approach to study interact- 
ing fermions in dimensions higher than one has been advanced by Castellani, Di Castro and 
Metzner |[24||. This method is based on Ward identities and leads to results that are very 



similar to those of Refs. 17 - 23 



Recently two of us have developed a functional integral method to bosonize interacting 
Fermi systems in arbitrary dimensions [25], which avoids the algebraic manipulations of 
commutators that is necessary in the operator approach. A similar method has been devel- 
oped independently by Frohlich and collaborators [p6| |, although in their work the physical 
content is somehow hidden in elegant mathematics. In the context of the one-dimensional 
Tomonaga-Luttinger model the functional bosonization technique has first been discussed 
by Fogedby [p7], and later by Lee and Chen [28]. Besides reproducing all of the results of 
Refs. jT7||- [[^J with a minimum of algebra JTB| , p5| , our functional bosonization approach has 
several other advantages. First of all, it can cope with retardation and non-locality in a very 
simple way, because path integrals are the natural language to describe these important 
many-body effects. Hence, the inclusion of electron-phonon interactions or the treatment 
of the retarded interaction mediated by transverse gauge-fields is straightforward in our ap- 
proach [|TIJ. A second important advantage of our approach is that it opens the way for a 
controlled calculation of corrections to the non-interacting boson approximation. As already 
mentioned, these non-Gaussian corrections exist in any dimension if the energy dispersion is 
not linearized at the Fermi surface. Some time ago Haldane pointed out in his seminal article 
on one-dimensional interacting Fermi systems [Q that the perturbation theory for the effec- 
tive interacting boson-problem, which arises from the bosonization mapping, is most likely 
well behaved. In this paper we shall verify this hypothesis not only in d — 1, but calculate 
explicitly the leading correction to the non-interacting boson approximation in arbitrary 
dimension. We shall show that this correction is closely related to the local field correction 
to the random-phase approximation (RPA) for the density- density correlation function. 

The rest of this paper is organized as follows: In Sec.|| we shall give a detailed description 
of our functional bosonization technique, and develop the general formalism that is necessary 
to calculate the corrections to the non-interacting boson approximation. Sec.|ITT| contains the 
calculation of the bosonized Hamiltonian and the density-density correlation function within 
the Gaussian approximation. The calculation of the leading correction to the Gaussian 
approximation will be presented in Sec.[TV]. We show that the non-Gaussian correction 
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terms to the bosonized Hamiltonian correspond to the local field corrections to the RPA, 
and calculate the hidden small parameter which determines the range of validity of the 
Gaussian approximation. Finally, in Sec.|V] we present our conclusions and discuss some 
open problems. 



II. GENERAL FORMALISM OF FUNCTIONAL BOSONIZATION 

In this section we shall describe our functional bosonization approach and prove that 
under certain conditions there exists for all d a remarkable cancellation between self-energy 
and vertex corrections, which implies that the Gaussian approximation becomes exact. 



A. The effective action for collective density fluctuations 

It is convenient to derive the bosonized Hamiltonian by first considering the imaginary- 
frequency density-density correlation function of our interacting many-body system, which 
is defined by 

II(q, iu m ) = 1%) = ^f Q dr J* dr'e-^-^ < T [p q (r)p_ q (r')] > , (2.1) 

where j3 = 1/T is the inverse temperature, p q = SkCk^k+q is the operator representing 
the Fourier transform of the total density, and T denotes time-ordering in imaginary time. 
Throughout this work q = [q, iu m ] is a collective label for wave- vector q and bosonic Mat- 
subara frequency u m = 2nmT. The average in Eq.[2.1| denotes thermal average with respect 



to all degrees of freedom in the system. An essential first step is the decomposition of Eq.[27T 



into contributions from the various patches. Using Eq. 1.5, we may write 



Tl( q )=J2n aa '(q) , (2.2) 

aa' 

where 

ir-'(g) = dr dr'e-^-r') < T [^( T )^' q ( T ')] > . ( 2 .3) 

We shall refer to 11(g) as the global density-density correlation function, and to T\. aa ' (q) as 
the patch density-density correlation function. IP a (g) can be represented as a functional 
integral over a Grassmann-field ip and an auxiliary field <fi" that mediates the interaction 
between the patch densities |T6|j25|j29|j3r 



W V J V {4)}V {cf) a } e -SiW-S2W,4> a }-s 3 {^} ' \ ■ ) 



where is now a quadratic functional of the Grassmann-fields, 

p« = J2<d a (k)4^ k+q . (2.5) 
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Here ipk is the Fourier component of the Grassmann-field, where k is again a collective 
label for wave-vector k and fermionic Matsubara frequency Cj n = 2ir(n + |)T. The three 
contributions to the action in Eq.E^l are 



SiM ^E4[-^+ai* , (2.6) 

s 2 {^,<p a } = *EE*v? > ( 2 - 7 ) 

q a 

s*m ^EEC'r^-/ , (2-8) 

Q act' 

where / is a dimensionless matrix in the patch labels, which is defined in terms of the 
coarse-grained Landau parameters f® a ' of Eq. |1.13| via 

^[L q r' = fT' ■ ( 2 - 9 ) 

The field 0" can be considered as the dual field of the composite fermionic field PRpB . 



We have intentionally chosen our notation for the actions Si, S 2 and S3 such that it matches 
the notation used by Feynman and Hibbs [31]: Si is the action for the matter degrees of 



freedom, S 2 describes the coupling between matter and a bosonic field that mediates 
the interaction between the matter, and S3 is the action of the bosonic field. In the case 
of the three dimensional Coulomb interaction the field 4> q can be identified with the scalar 
potential of electromagnetism in the Coulomb gauge, and it is straightforward to include 
also the transverse radiation field in our formalism JTH] . 

In order to introduce a collective bosonic field that represents the bosonized local density, 
we decouple S3 {4> a } by means of another Hubbard- Stratonowich transformation p^ , |25| , |2"8|| . 
The integration over the Grassmann-field ip is then quadratic and can be formally carried 
out, so that we can eliminate the composite fermionic field p® in favour of a collective bosonic 
field p~q. We arrive at the following exact expression for the patch density-density correlation 
function, 

8 fV\p a }e- s "^ a }p a J«' fDifje'^^^e-^W"} 

Yl aa (n) = — r-qrq J I r J ,^ 



where 



Q act' 

S kin {(f> a } =-Tr In [1-G V] . (2.12) 



Here the trace is over all wave-vectors and frequencies, and Go and V are infinite matrices 
in wave-vector and frequency space, with matrix elements given by 

[G ) w = S w G (k) , G (k) = - 1 - , (2.13) 
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(2.14) 



The subscript on Skin {4> a } indicates that this quantity is closely related to the bosonized 
kinetic energy. To be precise, below we shall show that the bosonized kinetic energy is given 
by the negative logarithm of the functional- Fourier transform of e 



s kin {<t> a } 



i.e. 



S kin {~p«} = -\n 



V{(f) a }e 



y , 



With this definition Eq. |2.10| can also be written as 

~ pC "' jO" 

rao' f„\ _ I" •> ~ *~r j - H-qHq 



U aa (q) 



/3 JV{p a }e- S -ff { P° !} n a n a ' 



leff 



{P a } 



with 



S eff {p a } = S int {p a } + S km {p a } 



(2.15) 



(2.16) 



(2.17) 



Note that the electron-electron interaction has been treated exactly, while the bosonized 
kinetic Skin{p a } can in general only be calculated approximately. The calculation consists 
of two steps. First of all, one should calculate the functional Skin {<p a }- I n praxis, this has 



to be done by expanding the logarithm in Eq. 2.12 and performing the traces in each term 
separately, 



Skin {0"} — £ Skin,n {0°} 
n=l 

S km ,n{<t> a } = -Tr \GoVT 

11 L J 



(2.18) 
(2.19) 



To obtain the bosonized kinetic energy Ski n {p a } it is necessary to perform the functional 
Fourier transformation in Eq. f2.15 . Of course, in general also this step has to be carried 



out perturbatively. If we ignore all terms with n > 3 in the expansion |2.18| , we obtain 



the Gaussian approximation. In Ref. [25] we have pointed out that in the one-dimensional 
Tomonaga-Luttinger model all higher order terms vanish identically, so that we have exactly 



— Tr In 



1-GoV 



Tr 



GnV 



+ -Tr 
2 



GnV 



[2.20) 



The functional Fourier transformation of Ski n {4> a } involves then a simple Gaussian integral, 
so that the bosonized kinetic energy can be calculated exactly. In the next subsection we 
shall show that in a certain physically relevant limit the remarkable cancellation between 
self-energy and vertex corrections that is responsible for the validity of Eq. [2.2(J| in d — 1 
happens also in higher dimensions. 



B. Generalized closed loop theorem 

Graphically, the trace in Eq. |2.19| can be represented as a closed fermion loop with n 
external Q -fields, see FigH Performing the trace in Eq.|2.19| we obtain 
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^n,n{0 a } = -£ £ U n (q iai ...q nan )<f>% •••</>% , (2.21) 

" gi...g n c«i. ..On 

where the dimensionless vertices U n are given by 

U n {q l a 1 ...q n a n )=5 qi+ ... +qnfi [^\ — £ £0 Qp i(k) 

\P/ n - P(l...n) k 

x e«^(k + qpj • • -e^(k + q Pl + . . . + qp n _J 

x G (k)G (k + g P J ■ • ■ G (k + q Pl + . . . + qp^) . (2.22) 

Here <5g 1+ ... +1?n> o denotes a Kronecker-5 in wave-vector and frequency space. We have used 
the invariance of Ski n ,n {0°} under relabeling of the fields to symmetrize the vertices U n with 
respect to the interchange of any two labels. The sum J2p(i... n ) * s over the n\ permutations 
of n integers, and Pi denotes the image of % under the permutation. Note that the vertices 
U n are uniquely determined by the energy dispersion ek — \i that defines the free fermionic 
action. The amazing point is now that there exists a physically interesting limit where all 
vertices U n with n > 3 vanish. This limit is characterized by the requirement that the 
following two approximations (Al) and (A2) become accurate: 



(Al): High density-limit or small momentum-transfer limit: Suppose that there exists a 
cutoff q c <C kp such that the contribution from fields 0° with |q| > q c to physical observ- 
ables becomes negligibly small. Because the fields <pg mediate the interaction between the 
fermions, this condition is equivalent with the requirement that the nature of the bare inter- 
action f ^ should be such that the resulting effective screened interaction (which takes into 
account the modification of the bare interaction between two particles due to the presence 
of all other particles) is negligibly small for |q| > q c . If this condition is satisfied, we may 
approximate in Eq. [2.22| 

Q a Pl (k)e «P 2 (k + qpj • • • Q ap - (k + q Pl + . . . + qp n _J 

w S ap i ap 2S ap i ap 3 ■ ■ ■ S ap i ap "Q ap i (k) , (2.23) 



because the k-sum in Eq. [2.22 is dominated by wave- vectors of the order of fc P . This ap- 



proximation is correct to leading order in q c /kp } and becomes exact in the limit q c /kp — > 0. 
Note that this limit is approached either at high densities, where kp — > oo at constant q c , 
or in the limit that the range q c of the effective interaction in momentum space approaches 
zero while kp is held constant. It follows that up to higher order corrections in q c /kp the 
vertex U n (qiai . . . q n a n ) is diagonal in all patch labels, 

U n (q iai . . . q n a n ) = 5 a ^ ■ ■ ■ 6 a ^U^ ( gi . . . q n ) , (2.24) 

with 

U:(q,...q n ) = 5 qi+ ... +qn J^] 1 £ £9 a (k) 

VP/ n - p(i...„) k 

x G {k)G {k + q Pl ) ■ ■ ■ G {k + q Pl + ... + ?Pn _J . (2.25) 
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Below we shall refer to the approximation |2.23| as the diagonal-patch approximation. It is 



important to note that at finite q c /kp this approximation can only become exact in d = 1, 
because in this case the Fermi surface consists of two widely separated points. Except 
for special cases |l3l , |l^ , |T6| , the covering of the Fermi surface in d > 1 involves at least 
some adjacent patches, so that in higher dimensions there exist always patches which can 
be connected by an arbitrarily small momentum transfer q. These " around-the-corner" 
processes are ignored within the diagonal-patch approximation (Al). 

(A2): Local linearization of the energy dispersion at the Fermi-surface: Suppose we 
linearize the energy dispersion locally by approximating « v a • q, see Eqs.|L9] and |1.10 . 



Shifting the summation wave-vector in Eq. [2.25| according to k = k a + q and defining the 
linearized Greens-function 

Gg(g)=Gg(q,»a w ) = _\ a - , (2.26) 

iLu n v • q 

Eq. |2.25| reduces to 



K{ qi ...q n )=5 qi+ ... +qn J 1 -) 1 £ E9 a (k° + q) 

\P J n. p(l ^ q 

x G a (q)G«(q + q Pl )---G a (q + q Pl + ... + q Pn _ 1 ) . (2.27) 



Having made the approximations (Al) and (A2), we are now ready to show that in 
arbitrary dimension the vertices U%(qi . . .q n ) with n > 3 vanish in the limit g c /A — > 0, 
so that in this limit the Gaussian approximation becomes exact! In the context of the 
Tomonaga-Luttinger model the vanishing of the U n , n > 3 has been called "weak closed 
loop theorem", and is discussed and proved in unpublished lecture notes by T. Bohr |32 



Under assumptions (Al) and (A2) the proof goes through in any dimension without changes. 
Note that the validity of (Al) and (A2) is implicitly built into the Tomonaga-Luttinger model 
by definition. The vanishing of U n for n > 3 is equivalent with the statement that the RPA 
for the density-density correlation function is exact in this model. This is due to a complete 
cancellation between vertex- and self-energy corrections. In d = 1 this cancellation has been 
discovered by Dzyaloshinskii and Larkin |33|]. Because the proof given in Ref. |33|] is rather 



cryptic and the lecture notes by T. Bohr [32] are not published, we briefly outline the basic 
features of the proof. It is perhaps simpler to formulate the proof in the space-time domain 



3~2| , |3"4|| , but for our purpose it is more convenient to work in momentum-space, because here 



the Fermi-surface and the patching construction are defined. The following two properties 
of our locally linearized Greens-function defined in Eq. |2.26| are essential, 



G a Q (~q) 



-G%(q) 



(2.28) 



G a Q ( q )G«(q + q') = G«(q') [G«(q) - G%(q + q')) 



(2.29) 



Note that Eq. |2.28| follows trivially from the definition, while Eq. |2.29| is nothing but the 
partial-fraction representation of the product of two rational functions. To show that the 
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odd vertices Uz,TJc,,... vanish, we only need Eq. |2.28| and the fact that the patch a in 
Eq. 2.27 has inversion symmetry with respect to k a , so that the domain for the q-sum is 
invariant under q — > — q. Then it is easy to see that the contribution from a given per- 
mutation (P1P2 ■ ■ ■ P n -iP n ) is exactly cancelled by the contribution from the permutation 
(P n P n _i . . . P2P1) in which the loop is traversed in the opposite direction. As already pointed 
out by T. Bohr [32|| , the vanishing of the odd vertices is a direct consequence of Furry's the- 
orem P5| . To show that the even vertices U n , n — 4, 6, . . . vanish, we use Eq. |2.29| n-times 
for the pairs 



GMG%(q + q Pl ) , 

G%{q + q Pl )G%{q + q Pl + qP2 ) , 

G%(q + q Pl + ... + q Pn _ 2 )G«(q + q Pl + . . . + q Pn _;, 
G%(q + q Pl + . . . + q Pn _ x )G%{q) , 

and take into account that we may replace q Pl + . . . + q Pn _ x = 
energy-momentum conservation. Then it is easy to show that 

G%{q)G%{q + q Pl ) ■ ■ ■ G a (q + q Pl + . . . + q Pn _ x ) 
= i [G%(q Pn ) - G«(q Pl )} G%(q + q Pl )--- G«(q + q Pl + 



[2.30) 



-q Pn because of overall 



+ 9p„-i) 



(2.31) 



Substituting Eq. |2.31| in Eq. |2.27| , shifting the summation label q — > q — q Pl + q Pn , and using 
the fact that we may rename q Pl <-> q Pn because we sum over all permutations, it is evident 
that the resulting expression vanishes. This argument is not valid for n = 2, because in 
this case Go^ft) — ^q(9pi) = 2GQ(gp 2 ) due to energy- momentum conservation. We shall 
discuss the vertex U 2 in detail in Sec JIII A| . Note that the shift q — > q — q Pl + q Pn affects 
also the patch-cutoff, O a (k a + q) — *• a (k Q + q — q Pl + q Pn ), but this leads to corrections 
of higher order in q c . Because we have already ignored higher order terms in q c by making 
the diagonal-patch approximation (Al), it is consistent to ignore this shift. 

The above proof of the generalized closed loop theorem uses the language of Greens- 
functions in the momentum-frequency domain. Such a formulation naturally emerges from 
our functional bosonization approach. After completion of this work we have learnt that 
independently from us W. Metzner has also realized that the closed loop theorem can be 
generalized to higher dimensions |3l| . His proof is based on the analysis of the commutation 
relations and equations of motion of the patch density operators defined in Eq. |l.ll . 



In fermionic language, the vanishing of the higher order vertices is due to a complete 
cancellation between self-energy and vertex corrections. This cancellation is automatically 
incorporated in our bosonic formulation via the symmetrization of the vertices U n . We 
would like to emphasize again that this remarkable cancellation happens not only in d — 1 
but in arbitrary dimensions. The existence of these cancellations in the perturbative 



calculation of the dielectric function of the homogeneous electron gas in d — 3 has already 
been noticed by Geldart and Taylor more than 20 years ago ||36|| , although the origin for 
this cancellation has not been identified. The generalized closed loop theorem discussed in 
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the present work gives a clear mathematical explanation for this cancellation to all orders 
in perturbation theory. 

It is important to stress that the cancellation does not depend on the nature of the 
external fields that enter the closed loop; in particular, it occurs also in models where the 
fermionic current density is coupled to transverse gauge-fields |[16|| . The one- loop corrections 
to the RPA for the gauge-invariant two-particle Greens-functions of electrons interacting 
with transverse gauge-fields have recently calculated by Kim et al. [071]. They found that 



at long wavelengths and low frequencies the leading self-energy and vertex corrections can- 
cel. In the light of the generalized closed loop theorem this cancellation is not surprising. 
However, the generalized closed loop theorem is a much stronger statement, because it im- 
plies a cancellation between the leading self-energy and vertex corrections to all orders in 
perturbation theory. 

A system where the approximations (Al) and (A2) are exact is the natural generalization 
of the Tomonaga model to higher dimensions ]T[|. If we follow Luttinger's construction |3j 
and remove the radial cutoff A by extending the locally linearized energy dispersion beyond 
the squat boxes (so that it is allowed to shift q) , we obtain an interacting higher dimensional 
Fermi system that is exactly solvable in precisely the same sense as the Tomonaga-Luttinger 
model in d — 1. As already mentioned, however, for finite q c /kp and realistic Fermi surfaces 
the condition (Al) can never be exactly satisfied in d > 1. 



III. GAUSSIAN APPROXIMATION 

In this section we shall calculate the density-density correlation function and the 
bosonized Hamiltonian within the Gaussian approximation. We also show that at long 
wavelengths the resulting bosonized Hamiltonian is identical with the Hamiltonian obtained 
from the operator approach |fL8| , p2]l . 



A. The first two vertices 

Within Gaussian approximation only the vertices U\ and JJ% are retained, and all higher 
order vertices are set equal to zero. From the previous subsection we know that the Gaussian 
approximation is justified in a parameter regime where (^41) and (A2) are accurate. Let us 
start with the vertex U\, which is given by 

^iM = ^^ 9 °( k ) - 1 c =<oN£ , (3.32) 

where 

iVi? = £e a (k)/(&) , (3.33) 

k 

is the number of occupied states in patch a in the non-interacting limit. Here /(e) = 
[e^ e + is the Fermi- function. Thus, 
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(3.34) 



Combining this term with the 0-p-coupling in the functional Fourier transform in Eq. |2.15 
we can write 



E E - s kin ,i m = * E E 



q a 



q a 



Thus, Skin,i {</> a } simply shifts the collective density field according to 

Pq Pg - <5,,0-/V . 



(3.35) 



(3.36) 



i.e. the uniform component is shifted. We conclude that the action Skin {p a } is actually a 
functional of the shifted field. For simplicity we shall from now on redefine the collective 
field according to Eq. [3.36| . Note also that the term q = in S int {p a } is usually excluded due 
to charge neutrality, so that the effective action S e ff{p a } can be considered as functional 
of the shifted field. 



The second-order vertex is given by 

1 



•n+*fl E [e ai (k)6 Q2 (k + qi )G (k)G (k + qi ) 



U 2 (qia 1 ,q 2 a 2 ) = -5, 

+ e a2 (k)Q a ^k + q 2 )G (k)G (k + q 2 )} . 
Performing the frequency sum we obtain 



U 2 (-qa,qa') = ^Uf(q) 



where 

nr'(g) 



2V 



E 



k L 



e a (k)e a '(k- q ; 



/(&-,) 



/(&) 



6c — £k-q — 



+ (a <-> a',q — > -g) 



Hence Sfci n ,2{0 a } is given by 



2/3 



9^9 



9 an' 



(3.37) 



(3.38) 



(3.39) 



(3.40) 



Note that Uq ' '(g) is the patch polarization in the absence of interactions, i.e. the exact 
result of Eq. |2.3| if we set / = 0. In this case our complicated transformations are of course 
not necessary and it is much simpler to calculate the density- density correlation function 
directly from its definition in terms of fermionic operators. For |q| <C kp the diagonal- 
patch approximation (Al) is justified, so that in Eq. |3.39j we may replace Q a {k)Q a (k - q) « 
S aa G a (k). To leading order in |q|/fcp we have therefore in any dimension 



U^\q)^6 aa 'U^q) , Tl«(q) = v c 



v Q • q — iu r 



(3.41) 
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where 



1 ON* 



V 



E© Q (k) 



(3.42) 



is the "local" density of states associated with patch a, and v a is the local Fermi velocity, 
see Eq. |1.10| . Note that Eq. |3.41| is valid for small |q|/fczr but for arbitrary frequencies. 



B. Density-density correlation function 



Within Gaussian approximation the action for the dual field 0" is quadratic, so that the 



functional Fourier transform in Eq. 2.15 can be calculated trivially, with the result 

4n{n = €,o+jEE r ^-/ > ( 3 - 43 ) 



q olol' 



where Sj^ n0 is a constant independent of the p"-field, and T" a ' is the propagator of the 



-field with respect to the Gaussian action Ski n ,2 {4> a }, 



'qYq> / Q - Oq+q'fi 1 q 



(3.44) 



Here < . . . >o denotes Gaussian average with probability measure e s ^,2{<t> a }^ j e f or an y 
functional F {4> a } 



(F{r}) 



JT>{(f) a } e- Sti ». 2 W"l 



(3.45) 



Note that T^ a ' is the inverse of the matrix U^—qa, qa') in the space spanned by the patch- 
indices, i.e. 



J2U 2 (-qa,qa')T a g ' a " = 5 C 



(3.46) 



Thus, Tq a ' is proportional to the inverse of the non-interacting patch polarization. Absorbing 

for simplicity the constant S^ n0 in the re-definition of the integration measure V {p a } and 
combining the contribution from the kinetic energy with the interaction contribution, we 
obtain for the bosonized action S e ff {p a } defined in Eq. [2.17 in Gaussian approximation 



Seff{p a } = ^T,[{lr'+^q a ' 



~P-q~Pq 



(3.47) 



q aa! 



Note that the inclusion of interactions is trivial, while the major difficulty lies in the repre- 
sentation of the kinetic energy as a functional of the collective density field. The Gaussian 
propagator of the p°-field is then 



fv{p a }e- s *ffV a yp«p«; 



V{p a }e~ s -ff^ 



Iq + r* 



(3.48) 
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where T_ q is a matrix in the patch-labels, with [I^] **' = T q a '. From Eq. |2.16| we obtain for 
the density-density correlation function 



U aa '(q) 



P 



V 



where H (g) is again a matrix in the patch-labels, with 



[n (g)r 



v L-« 



(3.49) 



(3.50) 



see Eqs. p.38| , p. 39 and |3.46| . The matrix elements of are the usual Landau-parameters 
defined in Eq. |1.13| , i.e. [/ ] aa ' = /° a . Eq. |3.49| is nothing but the RPA for the patch density- 
density correlation function. To obtain the standard RPA-result for the total density- density 



correlation function, we should according to Eq.|27| sum Eq. |3.49| over both patch-labels 



n(g) = Ef[no 1 (?)+i 



(3.51) 



For simplicity let us assume that [f 



f„ is independent of the patch-indices. Expanding 



[n l (g) + L q r l = n (g) - n (g)/ 9 n ( g ) + n (g)/ 9 n ( g )/ 9 n ( g ) 

and taking matrix elements, we see that Eqs. |3.51| and p.52 reduce to 

n (g) 



n( g ) 



i + / q n (g) 



where 



n (g) = £nr'(g) 



(3.52) 



(3.53) 



(3.54) 



is the total non-interacting polarization. We would like to emphasize that up to this point 
we have not linearized the energy dispersion, so that Eq. |3.53| is the exact RPA result for all 
wave-vectors, including the short wavelength regime. 



C. Bosonization of the Hamiltonian 



To make contact with the operator-formalism to bosonization |fL7|l - |23[| let us now derive a 
bosonic Hamiltonian that is in the limit of long wavelengths equivalent to our Gaussian action 
S e ff {p a } in Eq. |3.47| . The key observation is that the patch polarization is in the limit of high 
densities and long wavelengths (i.e. in the limit where the diagonal-patch approximation (Al) 
is correct) diagonal in the patch indices, and of the form given in Eq. |3.41| . From Eqs. |3.38| , 
3.39| and [3.46| it is then easy to see that the matrix elements of T q are given by 
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-pact 
3 



j (3 v" ■ q - iu^ 
v « . q 



|q| "C kF 



(3.55) 



Hence the Gaussian action in Eq. |3.47 can be written as 



Sejf in 



EE 



v" • q 



(3.56) 



Obviously the term proportional to iuj m defines the dynamics of the p a -field. We now recall 
that in the functional integral for canonically quantized bosons the coefficient of the term 
proportional to —iuj m should be precisely (3. Thus, to write our effective action in terms of 
a canonical boson field b^, we should rescale the p Q -field accordingly. This is achieved by 
substituting 



Pa 



(Vv a \v a ■ q|)3 [6(v Q • q)6° + 0(-v a ■ q)&t a 



(3.57) 



in Eq. p.56| . The G-functions are necessary to make the coefficient of —iuj m equal to (3 for all 
patches, because the sign of iui m in Eq. |3.56| depends on the sign of v" ■ q. Our final result 
for the bosonized action S {b a } = S e ff {p a {b a )} is 

S{b a } = /3EE ( va ■ q)(-^m)6f6° 



+ P[H km {b a } + H mt {b a }} 



(3.58) 



q a 



H*nt{b a } =^EE0K-q) V / |v Q -q||v Q '-q| 

9 aa' 



X 



©K' ■ q ) (j^'&X + 6; 
0(-v Q ' • q) (iT'h^l + r' a ^>; 



(3.59) 



(3.60) 



where /™' = \J v a v a> /™' are dimensionless couplings. The functional integral for the 6 a -field 

is now formally identical with a standard bosonic functional integral. The corresponding 

'rb 



inn where H kin and are 



second quantized bosonic Hamiltonian is therefore H b = H kin +H b 
simply obtained by replacing the bosonic fields b^ in Eqs. |3.59| , |3.60| by operators b^ satisfying 

[S°,6q,' f ] = 5 aa '5^. The resulting H b agrees with the bosonized Hamiltonian derived in 
Refs. [18|,|22| by means of an operator approach. 

Note, however, that the above identification with a canonical bosonic Hamiltonian is only 
possible in the limit of long wavelengths and high densities, so that our parametrization of 
the effective Gaussian action in Eq. [3.47| is more general. Moreover, for practical calculations 
the substitution in Eq. |3.57| is not very useful, because it maps the very simple form of 
S e ff{p a } given in Eq. |3.47| onto the complicated effective action S {b a } in Eqs. [3.58H3.6C 
without containing new information. 
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IV. BEYOND THE GAUSSIAN APPROXIMATION 



In this section we shall calculate the one-loop correction to the Gaussian approximation 
for the effective action of the collective density field. In this way we determine the hidden 
parameter which determines the range of validity of the Gaussian approximation. We also 
show that bosonization leads to a new systematic method to calculate corrections to the 
RPA for the dielectric function. 



A. General perturbative expansion of Sf-in {p a } 



In this subsection we develop a systematic perturbative method to calculate the bosonized 
kinetic energy Skin {p a } defined in Eq. |2.15| . This method is based on a loop- wise expansion 
of the functional Fourier transformation with the help of the linked cluster theorem. Defining 
S' kin {4> a } to be the sum of all non-Gaussian terms in the expansion of the effective action 
of the dual field (see Eq. j2.18p , 



sl m = e s k in,n m 

n=3 



we need to calculate 



e -Skin{p a } — g~ S £2l,0~' S: fe < si,2{£ a } /g-SLW"} 



where we have defined 



Q aa' 



and 



JV{r}F{r}e 



JV {(j) a } e^l" ^>- q Pg- S kin,2{<t> a } 

for any functional F {<ft a }. Performing the shift-transformation 



(4.1) 



(4.2) 



(4.3) 



(4.4) 



(4.5) 



it is easy to see that 

\ a' I 

In our case we have to calculate 



(4.6) 



(4.7) 
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Consider now the term of order in the expansion of S' kin {4> a } given in Eq. [4.1| . Clearly 

the substitution 0^ — > 4>g +iJ2 a ' ^q a Pq generates also a term of order (p a ) n which does not 
depend on the Q -field, and can be pulled out of the average in Eq.|4.7|. Let us denote this 
term by S { k °l n {p a }. From Eq.glH] it is easy to see that Sf2i n {P a } i s obtained by replacing 
0g ^«Ea' ^q a 'Pq in S kin)n {<j) a } , so that it is given by 

a' 

\ E E rM( qi a l ...q n a n )p%---p% , 

n qi...q n ai...a n 

where for n > 3 the vertices 1^°) are 

Tf{q l a l ...q n a n )=i n £ U n ( gi a[ . . . q n a' n )T a q '^ . . . T< a - , n > 3 . 



(4.8) 



(4.9) 



Obviously the Gaussian action S^!f n 2 {p a } in Eq.|4~3] is also of the form [4.8| , with 



r2 0) (»?iaig2a2) = S qi+q2j0 T 



92 



(4.10) 



The vertex U\ has been absorbed into the re-definition of p^ (see Eq. |3.36|) , so that 
S { k °Lm = 0. Defining 



Skin {P } — Skin,0 + E ^kin,n {P } > 
n=2 

s'L ir, P a } = s' kin {r + i E r"°y } - ^{i E r aQ y } 

o' a 1 

the general perturbative expansion for S k in {p a } is 



(4.11) 
(4.12) 



S km {p a } = S { k °l{p a }-\n 



i + E 



^=1 



A;! 



(4.13) 



According to the linked cluster theorem the logarithm eliminates all disconnected dia- 
grams, so that Eq. [4.13| can also be written as 



s kin {p a } = s { k i{p a }-i: 



k=l 



(4.14) 



where the subscript c means that all different connected diagrams should be retained |39 
From this expression it is easy to see that S k i n {p a } is in general of the following form 



(4.15) 



n=l 



where S k i n fi is a constant independent of the fields that cancels in the calculation of corre- 
lation functions, and for n > 1 
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S kin ,n{p a } = - E E T n {q 1 a l ...q n a n )p«l---P2 > ( 416 ) 

71 qi-..q n Ui...a„ 

where the vertices T n have an expansion of the form 

oo 

T n {q 1 a x ...q n a n ) = ^Y ( £\q l a 1 ...q n oi n ) . (4.17) 



k=0 



Here describes the interaction between n collective density fields that are generated 
from all diagrams in the linked cluster expansion |4.14| containing k internal loops of the 



-field. Note that the vertices in Eq. [4.9| do not contain any internal a -loops, and are 
therefore the tree-approximation for the exact vertices T n . Because each internal a -loop 
attached to a vertex U m reduces the number of external 0-fields by 2, it is clear that for k > 1 
the vertices can only be determined by vertices U m with m > n. Within the Gaussian 
approximation the vertices U m with m > 3 are set equal to zero, while the contribution 
from U\ can be absorbed into the re-definition of the q = 0-part of see Eq. |3.36| . Hence 
in Gaussian approximation 

T 2 (-qa,qa') « T^(-qa,qa') = Tf , (4.18) 



where T^ a ' is defined in Eq. |3.46| , and 



rW = , for n > 3 or k > 1 . (4.19) 

Although Ti = at the level of the Gaussian approximation, the higher order terms will in 
general lead to a finite value of Ti, which describes the fluctuations of the total number of 
occupied states in the patches. These terms do not contribute to correlation functions at 
finite q, but are certainly important for the calculation of the free energy. The formalism 
developed above is the starting point for a systematic calculation of corrections to the 
Gaussian approximation. 



B. Explicit calculation of the leading correction to the Gaussian approximation 

The leading correction to the Gaussian approximation is obtained from the one-loop 
correction in our effective bosonic theory, which amounts to a two-loop calculation at the 
fermionic level. Note that we have mapped the problem of calculating a two-particle Greens- 
function of the original fermionic model onto the problem of calculating a one-particle 
Greens-function of an effective bosonic model. The latter is conceptually simpler, because 
the symmetrized vertices U n and automatically contain the relevant self-energy and 
vertex corrections of the underlying fermionic problem. This will become evident below. 

At one-loop order, it is sufficient to truncate the expansion of the interaction part 
S' kin {(j) a } of the effective action of the Q -field in Eq.[0] at the fourth order, i.e. we may 
approximate 
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S'kinir} « S kin>3 {<j> a } + S fcm , 4 {0 a } 

= o E E ^3(giaig2tt2g3tt3)0g 1 1 05 2 2 C 3 

919293 «1«2"3 

+ 7 E E ^4(gia 1 g 2 a 2 g 3 a3g4a4)CC<aC > ( 4 - 20 ) 

91929394 «1«203"4 

where the vertices U 3 and C/4 are defined in Eq. [2.22| . According to the general formalism out- 
lined above, the bosonized kinetic energy Ski n {p a } is obtained by calculating the functional 
Fourier transform of Ski n {4> a }- Within one-loop approximation it is sufficient to retain only 
the term k = 1 in the linked cluster expansion Oj| , so that 

s ki n m = sit {n + (s'ur, p a }) , (4.21) 

where 

Skin {P a } = Skin,0 + g E E T P-qPq 

9 act' 

+ 7E E r£\ qi a iq2 a 2 q 3 a 3 )p%p%p% 

919293 «1«2«3 

+ 7 E E ^{qi^^a^a^lp^PZ , (4.22) 

91929394 01020304 

with 

rf ( qi a iq2 a 2 q 3 a 3 ) = -1 E t/ 3 (gia , 1 g2a 2 g3a 3 )r^ a T^ a2 rg°3 ; (4 . 23) 

a 'l a 2 a 3 

ri 0) (gi« 1 g 2 « 2 g 3 «3g4«4)= E f / 4(gi« / 1 g 2 « 2 g 3 « 3 g 4 al)r^ a T^ Q2 r"3-3 r <«4 (424) 



A a 2 a 3 a 'i 



The correction term due to one internal a -loop is 

(S'Li<P a , P a })o = SjSLja + S$ n>1 {p a } + >2 {p a } , (4.25) 



where 



5ffi,o = |E E ^(-^g^-gW^r^r^ 3 , (4.26) 

* gqr' 01020304 



= E r !V)Po , (4-27) 

o 



9 Qo' 



with 



ri 1} («)=zE E ^3(-?ai,ga 2 ,0a 3 )rj aoti rr' , (4.29) 
9 010203 

if^-?**, ?<*') = -3 E E ^4(-gai,?a 2 ,-g' a3 ,g'a4)Cr^'rp . (4.30) 

qi 01020314 
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Recall that the superscript (1) indicates that these terms contain one internal bosonic loop. 



Thus, within the one-loop approximation the constant in Eq. |4.15| is Ski n ,o = Sj^ n + 



and the vertices T n in Eq. [4.17| are approximated by 

r 1 (ga) = r! 1) (ga) , (4.31) 
T 2 (-qa,qa') = Tf +T^(-qa,qa') , (4.32) 

T 3 (q 1 a 1 q 2 a 2 q3<y3) = ^zHqi^i^^Qs^) , (4.33) 
r 4 (giaig2a2g3«3g4a4) = ri 0) (giaig 2 a2g3«3g4a4) , (4.34) 

and all T n with n > 5 are set equal to zero. The term with Ti can be ignored for a calculation 
of correlation functions at finite q, because it involves only the q = component of the density 
fields. Furthermore, at the level of a one-loop calculation we may also ignore the term Ig , 
because the Gaussian expectation value of a product of three /^-fields vanishes. Combining 
the relevant contributions from the kinetic energy with the interaction contribution, we 
finally arrive at the following effective action 

Seff {P a } = \ E E + V f] P-X + \ E E ifHza, qol)p%Pt 



2 _ 

9 aa' 9 aa' 

^7 E E r<?\q 1 a 1 q 3 a 2 q 3 a 3 q4a i )j%p%P%f% • (4.35) 

91929394 01020304 



We emphasize that this effective action is only good for the purpose of calculating the one- 
loop corrections to the Gaussian approximation. At two-loop order one should also retain 
the terms with T3 and Tq. The last two terms in Eq. [4.35| contain the one-loop corrections to 
the non-interacting boson approximation for the bosonized collective density fluctuations. In 
the limit of long wavelengths we may again write down an equivalent effective Hamiltonian 
of canonically quantized bosons by using the substitution |3.57| . However, we shall not even 
bother writing down this complicated expression, because this mapping is only valid at long 
wavelengths and high densities, and does not lead to any simplification. For all practical 
purposes the parametrization in terms of the collective density-field p a is superior. We shall 
now use this parametrization to calculate the leading correction to the bosonic propagator 
and determine in this way the hidden small parameter which controls the range of validity 
of the Gaussian approximation. 

C. First order self-energy and hidden small parameter 

Let us define a dimensionless proper (or irreducible) self-energy matrix E* via 



~o ~o 



PqPq' / - °q+q',0 



f„ + r q -r q 



(4.36) 



where the probability distribution for the average is determined by the exact effective action 
S e ff{p a }, see Eqs. [2.16| and |2.17| . From Eq. |3.48| it is clear that the self-energy S* contains 



by definition all corrections to the RPA. Introducing the exact proper polarization matrix 
Il*(q) via 
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[in?)]- 1 = [not?)]- 1 - ^s* , 

the exact total density-density correlation function can be written as 

n(?) = E[[[n*(?)r 1 + i 1 "' 



(4.37) 



(4.38) 



If all matrix elements of / are identical and equal to / q , we may repeat the manipulations 



in Eqs. p.51| - |3.53| , so that Eq. |4.38| reduces to |3 

n*(g) 



11(g) 



i + / q n*(g) ' 



where 



n*(g) = £[rr(g)p' 



(4.39) 



(4.40) 



is the total proper polarization. With the help of the dielectric function 

e(g) = l + / q lT(g) , 

Eq. |4.39| can also be written as 

n*( g ) 



n( g ) 



(4.41) 



(4.42) 



To first order in an expansion in the number of bosonic loops, we simply have to add the 
two diagrams shown in Fig.|j. Because we have symmetrized the vertices, the loop diagram 
in Fig.f| has a combinatorial factor of three, so that to first order E* ~ TJ^\ with 



-T^i-qa^qa') — r 4 0) (-?«, -q'a 3 , q'a 4 ) 



(4.43) 



Using the definitions of r^ 1 "* and T^' (see Eqs. [4.30| and 4.24 ), it is easy to show that Eq. |4.43 
can also be written as 



[£«]<-' = 3£ E U 4 (-qa l ,qa 2 ,-q'a 3 ,q'a 4 )r^T^ a ' 



n' Q1Q20304 



X 



Eg/ — T q > 



-1 



«4«3 



A simple manipulation of the matrix in the last line gives 

-i_ ~. r _,, r i 



r - r 
— g —1 



r_ rRPA 



(4.44) 



(4.45) 
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where the matrix f^ PA is defined by 



RPA 



/Ji + n (g)^ 



-1 



It follows that Eq.[4.44| reduces to 



[£ (i) r ' = 3 ^ £ E T^T^'U 4 (-qa l} qa 2} -q'a 3} q'a,) [Q 



(RPA 

-g' 



O4O3 



(4.46) 



(4.47) 



Note that is proportional to the RPA-screened interaction and vanishes in the non- 
interacting limit, as it should. Eq. 4.47] is the general result for the leading correction to the 
Gaussian approximation in an arbitrary patch-geometry, and for an arbitrary interaction 

matrix f . 
—9 

To make further progress, we shall now restrict ourselves to the high-density limit and 
assume that the interaction is negligibly small for wave-vectors larger than a cutoff q c <C 
min{A, A} <C hp. Choosing also the magnitude of the external wave-vector q in Eq. |4.47 
small compared with the characteristic size of the patches, the diagonal-patch approximation 
(Al) is justified, so that £ and f/ 4 (— ga 2 , —q'a 3 , q'a^) are diagonal in the patch-indices, 
see Eqs. p.55| and |2.24j . Then Eq j4.47| reduces to 



Miliar' _ @ / v a • q - iiO m \ 

[ ^ J ~ 6 v« ■ q ) 9 

where the dimensionless function A°t is given by 

3/3 



(4.48) 



A a _ 

q Vv a 



T,Un-q,q,-q',q') f/ A 



(4.49) 



We conclude that in the high-density- and long wavelength limit we have to first order in 
the screened interaction 



■paa 



a a wa' 



(4.50) 



with 



1 - A«)v<* • q - (1 - 2A«)iu m - A«^£ 



q v Q q 



U* a (q) 



v a v a • q 



(4.51) 



Comparing this expression with Eq. |3.55| , we conclude that the non-interacting boson ap- 
proximation is quantitatively correct provided the condition 



\A«\ « 1 



(4.52) 



is satisfied for all patches a, because then the corrections to the propagator of the locally 
defined collective density field are small. 
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For simplicity, let us now assume that all matrix elements of / are identical. Using 



—q 



the same procedure as in Eq. |3.52| , it is then easy to show that [f RPA } aa = f RPA is also 



independent of the patch label. From the general definition of the vertices U n in Eq. |2.22| we 
obtain in this case 

A* = E ©°(k) {G (k)rtp(k)G (k)[G (k + q) + G (k - q)} 

+ \G (k)[T^(k, q)G Q (k + q) + T?(k, -q)G {k -</)]} , (4.53) 

with 

^F\k) = ^yEf q r A G (k + q') , (4.54) 
T { p(k,q) = ^-Y,fT AG o(k + q')Go(k + q' + q) . (4.55) 

Note that A°^ q = due to the symmetrization of the vertex U4 with respect to the in- 
terchange of any two labels. It is now obvious that the vertices of our effective bosonic 
action automatically contain the relevant self-energy and vertex corrections of the under- 
lying fermionic problem [36|. The first term in Eq. [4.53| corresponds to the two self-energy 



corrections to the non- interacting polarization bubble shown in Fig.|5] (a) and (b), while the 
last term is due to the vertex correction shown in Fig.[5] (c). 

In order to determine the range of validity of the non-interacting boson approximation, 
we have to calculate the dependence of A® on the various parameters in the problem. In 
the limit of long wavelengths and low energies, it is to leading order in |v Q • q| and \u m \ 
consistent to replace in Eq. |4.51| A^ — > Aq. Actually, the q — > limit of should be 



taken in such a way that the ratio is held constant, because in this case we obtain the 
low-energy behavior of A^ close to the poles of the Gaussian propagator. However, since we 
are only interested in the order of magnitude of A^ for small \u m \ and |q|, it is sufficient 
to consider the "q-limit" Aq = ]im n ^ o \lwi u , m _> A"]. For simplicity, let us now ignore the 
frequency-dependence of the RPA-interaction. This amounts to the static approximation 
for the dielectric constant, which seems reasonable to obtain the correct order of magnitude 
of Aq. The "q-limit" is obtained by setting q = under the summation sign and performing 
the Matsubara sums before doing the wave-vector integrations. For /3->oowe obtain 

^=^Ee»(k)/r{/& + -,-)^/&) + |/& + ,.)|/&)} • (4-56) 

Because the k-sum extends over the entire box K a and by assumption the q'-sum is cut 
off by the interaction at q c <C min{A, A}, we may set £k+q' ~ 6c in the Fermi-functions of 
Eq . |4.56| . Then the summations factorize, and we obtain 



4" 
^0 



J_ \ ' fRPA 
v q 



r d^(o mono + norm , (^v 
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where v a (£,) is the energy dependent patch density of states 

1 

k 

Note that from the definition of v a in Eq. |3.42 



^(0 = y£©°W(e-6O • ( 4 - 58 ) 



/oo 
dfr a (0[-f(0) . (4.59) 
-oo 

Integrating in Eq. |4.57| by parts, we obtain in the limit (3 — > oo 

/oo 
^(0 [/(0/"(0 + /'(o/'(0] = (4-60) 
-oo 

r°° d 1 <9v a 



Because by assumption / PPA becomes negligibly small for |q| > q c , the first factor in Eq. |4.57 
is for V — ► oo 

^E/r^/o flP V • (4-62) 
q 

Ignoring a numerical factor of the order of unity, the final result for Aq can be written as 

d fRPA 

v4«~Mo c a , (4.63) 



where the dimensionless parameter 



_d 2 NJ_ 



v a dpi 



is a measure for the /oca/ curvature of the Fermi surface in patch K%. For sufficiently 
smooth Fermi surfaces the patch density of states v a is proportional to A 1 . However, 
the cutoff dependence cancels in Eq.Ojl, because it appears in the enumerator as well as 



in the denominator. Therefore C a is a cutoff-independent quantity. If we linearize the 
energy dispersion in patch K", then the local density of states is replaced by a constant 
independent of fi, so that C a vanishes in this case. Then there is no correction to the 
Gaussian approximation. Of course, we already know from the closed loop theorem that 
the Gaussian approximation is exact if in addition to the diagona-patch approximation (Al) 
the energy dispersion is linearized for all patches. As usual, we introduce the dimensionless 
interaction F^ PA = ufQ PA 1 which measures the strength of the potential energy relative 
to the kinetic energy. Note that for regular interactions F^ PA m y+^- Using the fact to 
the global density of states is in d dimensions proportional to kfr 2 , we conclude that in 
the high-density limit the Gaussian approximation is quantitatively correct provided the 
condition 
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\A$\ = (J^j \F RPA \\C a \ < 1 (4.65) 

is satisfied. This is the main result of this section. 

The appearance of three parameters that control the accuracy of the Gaussian approxi- 
mation has a very simple intuitive interpretation. First of all, if at all points on the Fermi 
surface the curvature is intrinsically small (i.e. \C a \ <C 1 for all patches a) then the cor- 
rections to the linearization of the energy dispersion are small, and hence the Gaussian 
approximation becomes accurate. Note that in the one- dimensional Tomonaga-Luttinger 
model C a = 0, because the energy dispersion is linear by definition. However, for realistic 
energy dispersions of the form = ^— the dimensionless curvature parameter C a is of the 
order of unity. But even then the Gaussian approximation is accurate, provided the nature 
of the interaction is such that it involves only small momentum transfers. This is also intu- 
itively obvious, because in this case the scattering processes probe only a thin shell around 
the Fermi surface and do not feel the deviations from linearity. Finally, it is clear that also 
the strength of the effective interaction should determine the range of validity of Gaussian 
approximation, because in the limit that the strength of the interaction approaches zero all 
corrections to the Gaussian approximation vanish. 



Finally, we would like to discuss a subtlety associated with the choice of the radial patch 
cutoff A and the picture of having integrated out all degrees of freedom outside a thin 
shell around the Fermi surface. Consider the three dimensional Coulomb interaction in the 
homogeneous electron gas, which in the static screening approximation corresponds in the 
regime |q| <C kp to 

K PA = ^—, , (4-66) 
M q -|- k 

where k = [47re 2 z/] 1//2 is the Thomas-Fermi screening wave-vector. For |q| « kp the Coulomb 
potential reduces to a constant of the order of fr- If we introduce a radial patch cutoff A ~ « 
and define our bare model such that it contains only the degrees of freedom close to the Fermi 
surface (see Sec.|), then we see from Eq. [4.66| that we must choose q c ~ k and F PPA « 1, so 
that Eq. [4.65| gives 

AJwfAV , A«k , (4.67) 



k 



F 



where we have ignored numerical constants of the order of unity. Note that C a = 0(1) for 
free electrons in d — 3. On the other hand, if we do not introduce a radial patch cutoff 
but treat also the high-energy degrees of freedom explicitly, then our simple estimate in 
Eq. |4.65| is not valid, because F PPA in Eq. [4.66| does not fall off sufficiently fast, so that the 



summations in Eq. f4.56| do not factorize. Substituting Eq. |4.66| into |4.56| , it is not difficult 
to show that in this case the value of Aq is essentially determined by the short-wavelength 
regime k < |q| < kp, and that it is given by 

~ (t~) 2 , A = oo . (4.68) 



F 
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Although in both cases the corrections to the Gaussian approximation are small at high 
densities (where k -C kp), the corrections are smaller by a factor of n/kp in the model with a 
radial patch cutoff. This is evidently due to the fact that in the cutoff model the high-energy 
degrees of freedom have already been taken into account implicitly via the re-definition of 
the parameters in the Gaussian approximation, so that only the much smaller number of 
degrees of freedom in a thin shell around the Fermi surface can give rise to corrections. 
However, our approach does not depend on the linearization of the energy dispersion, so 
that the radial cutoff can be removed and also the high-energy degrees of freedom can be 
explicitly taken into account. This enables us to use our bosonization approach as a basis 
to develop a new systematic method of calculating corrections to the RPA. We shall discuss 
this method in the following subsection. 



D. Calculating corrections to the RPA via bosonization 



The corrections to the RPA are usually expressed in terms of the dynamic local field 
factor g(q), which is defined by writing the exact irreducible polarization (see Eq. |4.4U| ) in 
the form 



l-H 



n*(g) 



n (g) 



l-g(q)U (q) 
Note that this equation can also be written as 

miT 1 = WIT 1 - 9{q) 



(4.69) 



(4.70) 



which has the structure G^ 1 = Gq 1 — E, i.e. it resembles the Dyson equation for the single- 
particle Greens-function of a bosonic problem, with the proper polarization and the local 
field factor playing the role of the exact Greens-function and the irreducible self-energy. 
This analogy is well known PS]- although it seems that it has not been thoroughly 
exploited as a guide to develop systematic methods to calculate corrections to the RPA. 
Such a method naturally emerges from our bosonization approach. Defining a local-field 
correction matrix 



our matrix Dyson equation 4.37 for the patch proper polarization takes the form 



(4.71) 



(4.72) 



Comparison with Eq. [4.70| shows that the matrix elements [£*] QQ of the irreducible self- 



energy of our effective bosonic model can be identified physically with generalized local field 
corrections [g{q)] aa , which differentiate between the contributions from the various patches. 
For simplicity let us assume that the diagonal-patch approximation (Al) is justified, so that 
Eq. |4.72| reduces to an equation for the diagonal elements, which can be written as 
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U* a (q) 



ng(g) 



i-s°(?)no tt (?) 

^U^q)+U^ q )g a ( q )U^q) + ... . (4.73) 

Here Tl* a (q) = [T\*(q)] aa , g a (q) = [g(q)] aa , and IIq (q) is at long wavelengths given in Eq. |3.41| . 
The difference between our approach, which is based on the perturbative calculation of 
the inverse proper polarization, and the naive perturbative approach is now clear. In the 
latter method the corrections to the proper polarization bubble are determined by direct 
expansion of II* (q) in powers of the interaction [[5]]. Such a procedure does not correspond 



to the perturbative calculation of the self-energy, but is equivalent to a direct expansion 
of the Greens-function. To first order, only the leading correction in the expansion of the 
Dyson equation (the second line in Eq. |4.73|) is kept in this method, so that the total proper 
polarization is approximated by 



H*(g) « £ 



v v • q 



v a ■ q 



lOJr. 



+ v a A a q 



Ml) + E {G (k)^\k)G (k)[G (k + q) + G (k - q)} 



1 



G {k)[T<p{k, q)G (k + q) + T^(k, -q)G (k - q)] 



(4.74) 



see Eqs. [4.53[ - [4.55| . In the three dimensional Coulomb-problem, the correction term in Eq. |4.74j 
has been discussed in Refs. [36,40]. Note, however, that these authors evaluate the fermionic 
self-energy \k) and vertex correction Tp(k, q) with the bare Coulomb interaction. Holas 



et al. [fH]] find that the expansion in Eq. |4.74| leads to unphysical singularities in the dielectric 
function close to the plasmon-poles. 

In contrast to the direct expansion of II* (q) in powers of the bare interaction, in our 
method we first calculate the irreducible self-energy of the effective bosonized Hamiltonian 
in powers of the RPA-screened interaction, and then re-sum the perturbation series via the 
Dyson equation. The crucial point is that the problem of calculating corrections to the RPA 
can be completely mapped onto an effective bosonic problem: our functional bosonization 
method allows us to explicitly construct the interacting bosonic Hamiltonian. Once we 
accept the validity of this mapping, we can use standard many-body theory for bosonic 
systems, which leaves us no choice: The corrections to the propagator of this effective 
bosonic theory should be calculated by expanding its irreducible self-energy E* (i.e. the 
inverse Greens-function) in the number of internal bosonic loops, and then solving the Dyson 
equation. A similar re-summation has been suggested in Refs. [ 40, 42, 43 1, but it is not so 



easy to justify this procedure at the fermionic level. Our bosonization approach provides 
the natural justification for this re-summation. The unphysical singularities [[4(]] that are 
encountered in the direct perturbative approach might be understandable from the point of 
view of bosonization: they are artificially generated because one attempts to calculate the 
bosonic single-particle Greens-function by direct expansion. This expansion is bound to fail 
close to the poles of the Greens-function, i.e. close to the plasmon poles. 



Based on the insights gained from our bosonization approach, we would like to suggest 
that corrections to the RPA should be calculated by expanding the generalized local field 
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corrections g(q) in powers of the RPA-interaction and then solving the bosonic Dyson equa- 



tion [4.72| . We suspect that in this way unphysical singularities in the dielectric function can 
be avoided. In the high-density limit we obtain from the first line in Eq. f4.73| for the total 
proper polarization 



n*(g) = £ 



l-A; 



:V a q 



v Q -q-^ m (l-^ 



_ (iu) m ) 2 r> Q 
v Q q q 



(4.75) 



with 



Bt 



Act 



(4.76) 



A detailed analysis of Eq. |4.75| and the resulting dielectric function e(q) = 1 + / q II*(g) in 
various dimensions requires a careful analysis of the analytic properties of the function 
and extensive numerical work, which will be presented elsewhere. Here, we would like to 
restrict ourselves to the discussion of the compressibility K p , which can be obtained as the 
"q-limit" of the exact proper polarization (compressibility sum rule) ||39|| , 



lim 

q^O 



lim Il*(q, iu % 



(4.77) 



As already mentioned, for sufficiently smooth Fermi surfaces the dependence of A^ on the 



patch- index a can be ignored. Setting A = limq^oflim^^o A"], we obtain from Eq. |4.75 



K n 



v 



1-A 

If we combine this expression with the estimate for A given in Eq. |4.68| , we obtain 

v 



) 2 ' 



(4.78) 



(4.79) 



where 7 is a numerical constant of the order of unity. Eq. |4.78| agrees with the expression 
given by Singwi et al. [44]], who present also numerical results for 7 as function of the 



density of the electron gas. Note that according to Eq. [4.78| the compressibility diverges and 
the system becomes unstable for Aq — > 1. Eq. [4.78| should be compared with the expression 
for the compressibility resulting from approximation |4.74| , K p « v{\ + A ). This is the 



leading term in the expansion of Eq.4.78 and does not predict any instabilities 



V. CONCLUSIONS 

The main results of this work can be summarized as follows: 

(1) We have developed a functional integral method which allows us to bosonize in- 
terating fermions in arbitrary dimensions. In general, bosonization maps a system of in- 
teracting fermions onto an effective system of interacting bosons. The corrections to the 
non-interacting boson approximation can be calculated in a systematic way. 
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(2) There exists a physically interesting limit where the Gaussian approximation for the 
effective bosonic Hamiltonian becomes exact. This limit is characterized by the requirements 
that the diagonal-patch approximation (Al) and the local linearization (A2) are correct. 
If the approximations (Al) and (A2) are made and the radial cutoff A is removed, the 
resulting interacting Fermi system is exactly solvable in any dimension due to an exact 
cancellation of self-energy and vertex corrections (generalized closed loop theorem). In the 
one-dimensional Tomonaga-Luttinger model the condition (A2) is satisfied by definition, 
while (Al) is satisfied as long as the two Fermi points are sufficiently separated, so that 
they cannot be connected by the maximal momentum transfer of the interaction. Because 
for realistic Fermi surfaces in d > 1 there exist always neighboring patches that can be 
connected by arbitrarily small momentum transfers, (Al) is only approximately correct in 
higher dimensions. 

(3) If the diagonal-patch approximation (Al) is made but the energy dispersion is not lin- 
earized, we have explicitly calculated the leading correction to the Gaussian approximation, 
and have shown that the effect of interactions between the bosons is negligible as long as 

the parameter f|M l-^cf^l is small compared with unity. The dimensionless quantity 
C a is proportional to the derivative of the density of states with respect to the chemical 
potential, and vanishes if the energy dispersion is linearized. The origin of an additional 
small parameter C a in the correction to the RPA lies in the cancellation between self-energy 
and vertex corrections, as described by the generalized closed loop theorem. 

(4) Our bosonization approach maps the calculation of the density-density correlation 
function of the Fermi system on the problem of calculating the single-particle Greens- 
function of an effective bosonic model. This mapping explains the origin of unphysical 
singularities that arise in direct perturbative expansions of the proper polarization, and 
strongly suggests that corrections to the RPA should be calculated by expanding the gen- 
eralized local field corrections in powers of the RPA-interaction. 



The present work is the first step to go beyond the non-interacting boson approximation. 
We have calculated the bosonized Hamiltonian as well as the density-density correlation 
function beyond the Gaussian approximation. The problem of calculating the non-Gaussian 
corrections to the leading bosonization expression for the single-particle Greens-function 
remains open. This calculation is more difficult, because it involves the solution of a partial 
differential equation [fT6|, p5| , p8| ,fi5|| . If the energy dispersion is linearized, £° ~ v a • q + 0(q 2 ^ 



this differential equation is linear and first order, and can be solved exactly. Inclusion of non- 
linear terms of order q 2 and higher in the expansion of £° lead to higher spatial derivatives, 
which are not so easy to handle. Following the method used by Khveshchenko and Stamp 
H3fl , it might be possible to set up some kind of perturbation theory for the solution of this 



higher order partial differential equation, but so far we have not pursued this possibility. 

It is tempting to take the form of the Greens-function within the Gaussian approximation 
as a basis to speculate on the effect of the non-Gaussian terms. Within the Gaussian ap- 
proximation the general bosonization result for the real-space imaginary-time single-particle 



Greens-function associated with patch a can be written as p5] 



G a (r,r) = G£(r,r)e QQ(r ' r) , (5.1) 
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where Gg(r, r) is the non-interacting Greens-function, and 

Q a (r 1 r) = R a -S a (r,r) , R a = lim S a (v,r) , (5.2) 

r,r— »0 

1 / cos(q • r - uj m T) 

s "^-w^ - • (5 - 3) 

Note that this expression depends exclusively on the RPA-screened interaction matrix f^' PA , 
as defined in Eq. [4.46| . Due to the linearization of the energy dispersion the non-interacting 



real-space Greens-function Gq(v,t) is proportional to 5^ 1 \r±), where r_|_ is the d — 1- 
dimensional projection of r perpendicular to v a . Therefore at the level of the Gaussian 
approximation we may replace 

Q a (v,r) ^Q a (r^ a ,r) (5.4) 
in Eq.|5.1|, where v a is a unit vector in the direction of v Q and r» = r ■ v a . If the non-linear 



terms in the energy dispersion are taken into account, we expect that the RPA-interaction is 
replaced by the exact effective interaction, which takes the local field corrections to the bare 
polarization bubbles into account. Furthermore, the term v Q • q in the denominator of Eq.|5.3| 



should be replaced by the full dispersion and the substitution |5J] cannot be made any 
more. However, non-Gaussian terms will certainly also give rise to further modifications, so 
that it is not clear if the Greens-function retains at least approximately the above simple 
structure. 

Finally, it would be interesting to calculate the corrections to the RPA by means of the 
method outlined in Sec. [lVDl . This method is different from a straight-forward expansion of 



the proper polarization, because it is based on the loop-wise calculation of generalized local 
field corrections associated with the patches. The calculation of corrections to the RPA has 
been discussed intensely in the 70's, but today the interest in this problem has waned. Our 
bosonization approach sheds new light on the various approximations that can be found in 
the literature, and perhaps will revive the interest in this problem. 
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FIGURES 



FIG. 1. Graph of a squat box Kj[ A with patch cutoff A and radial cutoff A in three dimensions. 
The vector k a points on the Fermi surface to the center of patch K%. 

FIG. 2. The boxes for a spherical Fermi surface in d = 2. The dashed lines mark the 
corresponding squat boxes A that contain only the degrees of freedom in the vicinity of the 
Fermi surface. 



FIG. 3. Feynman-diagram representing Tr[GoV] n , see Eqs. 2.19| and 2.21. The lines with arrows 
denote fermionic Greens-functions, while the wavy lines denote external ^"-fields. 



FIG. 4. Leading contributions to the bosonic self-energy, see Eq.4.43. The dark circle denotes 
; vertex defined in Eq. 4.3C and the shaded square represents the vertex defined in 
The dashed arrows denote the collective density field p a , and the dashed loop is the 



Eq.4.24 



Gaussian propagator of the /5 a -field, see Eq. 3.48 



FIG. 5. Leading local-field corrections to the polarization bubble, (a) und (b) are the leading 
self-energy corrections (see Eq.4.54), while (c) is the leading vertex correction (see Eq j4.55 ). The 
thick wavy line denotes the RPA-interaction, as defined in (d). The thin wavy line in (d) represents 
the bare interaction. 
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